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Abstract 

We study boundary effects in a linear wave equation with Dirichlet type conditions in a 
weakly curved pipe. The coordinates in our pipe are prescribed by a given small curvature with 
finite range, while the pipe's cross section being circular. Based on the straight pipe case a 
perturbativc analysis by which the boundary value conditions are exactly satisfied is employed. 
As such an analysis we decompose the wave equation into a set of ordinary differential equations 
perturbatively. We show the conditions when secular terms due to the curbed boundary appear 
in the naive peturbative analysis. In eliminating such a secularity with a singular perturbation 
method, we derive amplitude equations and show that the eigenfrequencies in time are shifted 
due to the curved boundary. 

PACS: 02.60.Lj, 02.30.Mv, 02.40.Hw 

1 Introduction 

In a spatially extended system a perfunctory theoretical description to explain experimentally 
observed phenomena could be done by ignoring or simplifying boundary conditions. However 
in reality, the observed phenomena in finite domain are generally affected by spatial boundary 
conditions and are non-trivial. A reason why boundary conditions are ignored or simplified is the 
lack of mathematical tools. It has been difficult to study such boundary effects in a systematic 
manner. We then need to develop a theoretical framework that enables us to describe some effects 
of spatial boundary conditions. A first attempt might rely on a perturbative analysis, and the 
small parameter is related to the magnitude of deformation from trivial spatial boundary. In 
dynamical systems theory there are a variety of useful methodologies for dealing with perturbed 
systems pp. With those methodologies differential geometry provides powerful mathematical tools 
for it [2]. Beyond some existing works on this context, such as electromagnetic waves in a curved 
pipe[H H], quantum eigenstates of a curved nanowire[5], one would like to know the higher order 
corrections in the magnitude of deformation of boundaries. To clearly see what can be observed 
including higher orders due to non-trivial spatial boundary conditions one needs to have a simple 
model. As such a model, we consider a linear wave equation which is widely studied in physical 
sciences, for example, fluid dynamics, electromagnetism and high energy physics. In this paper 
the linear wave equation, or strictly speaking the classical complex Klein-Gordon equation, with 
Dirichlet type conditions is studied based on the language of differential geometry with the use of 
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perturbation methods. There we will derive amplitude equations. In fact the idea of amplitude 
equations is often used in order to study weakly perturbed systems in nonlinear science and give 
benefits. We will see that even in the linear equation a mode-coupling phenomenon occurs due 
to the weakly deformed boundary, and then perturbative correction terms for the mode-amplitude 
will be obtained using a singular perturbation method. 

2 Coordinate system, co-frame and Laplacians 

First we give the expression for the coordinate system adapted to our curved pipe. The z-coordinate 
for our curved pipe which we consider is assumed to be planer and prescribed by the curvature 
which is small. Accordingly we do not consider the spatial curve, and our curve is nearly straight. 
This assumption on curvature makes the perturbative analysis effective. The two-dimensional cross 
section at any z in our pipe is assumed to be circle whose radius is a ( See FigQ]). 
We denote the given curvature by 

k(z) = eK (z), (1) 

where |/«(z)a| <C 1 and e is the small parameter. 

To describe the geometry of boundary we use the following metric tensor 

g = e 1 <g> e 1 + e 2 <g> e 2 + e 3 ® e 3 , 

where {e 1 , e 2 , e 3 } is the co-frame to be specified later, and our space will be denoted by U. The 
canonical volume elements are defined as 

#l = e 1 Ae 2 Ae 3 , #1 = e 1 A e 2 . 

Then the Hodge maps # are defined : 

# : YA q U — > rA 3 " 9 U, # : TA q U -> TA 2 ~ q U, 

acting on { e , e 2 , e 3 } and { e , e }, respectively. Here TA q U denotes the set of (/-form fields on U. 

For our weakly curved pipe the adapted co-frame is derived from the use of Frenet frame. The 
space-curve is given in terms of Euclidean position vector C{z) then the points in the interior can 
be written as 

r(z, r, 6) = C(z) + xi(r, 6) n(z) + x 2 (r, 6) b(z). 




z=0 z=L 



Figure 1: Curved pipe 
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Here x\ and X2 can be chosen as x\ = rcos#,;c2 = rsin# with < r < a and < 9 < 2tt. In 
addition, n is the Frenet normal and b is the Frenet bi-normal vectors to the curve, the tangent 
vector t is given by dC /dz. The relations between them are describes as the Frenet-Serret formulae: 

d d d d 

— C = t, —t = kti, — n = —ki + to, — o = — rn, 

dz dz dz dz 

with t(z) being the Frenet torsion, which is zero in our case. Since the infinitesimal deviation of r 
is written with 6x1,6x2 and 6z that are infinitesimal deviations of xi,x% and z, 

6r = n{5x\ — x^jbz) + b(6x2 + x\t5z) + i(l — kx±)6z, 

one can choose a convenient orthonormal co-frame with r = for the interior domain 

{ e 1 = dr e 2 =rd9, e 3 = (1 — €Kq(z)t cos 9)dz }, 

Taking the limit e — > above one has the adapted co-frame for the straight pipe 

{ e 1 = dr e 2 = rd9, e 3 = dz }. 

These are the usual cylindrical ones. 

2.1 Laplacians 

The Laplacian on U for differentiable functions is 

#d#d, (2) 

where d : TNPIA — > FA P+1 U is the exterior differentiation in a patch with coordinates in IA. So, for 
any function of t,z,r,9 one has 

df(t, z, r, 9) = ^fdr + %d9 + jj-dz. 

or ro9 az 

The explicit form of ([2]) for functions is calculated to be 

,,,,,,( 1 d 2 K,'(z)r cos 9 d 

#d#d=i j- -——-— — + 



(1 — n{z)r cos 9) 2 dz 2 (1 — n(z)r cos 9) 3 dz 
d 2 d k(z) cos 9 (9 1 ( d 2 K(z)sin9 d 



dr 2 rdr 1 — n(z)r cos 9 dr J \r 2 d9 2 r(l — n(z)r cos 9) 89 

where ' denotes the differentiation with respect to z. When k = this expression of the Laplacian 
reduces to 

d 2 d 2 d d 2 

• • -T • "7777777 • (3) 



r=0 dz 2 dr 2 rdr r 2 d9 2 



A(°) -=#d#d 

Higher order corrections are obtained by substituting k(z) = €Kq(z), (H|), and expanding it in e as 

#d#d = A (0) + eA (1) + e 2 A (2) + • • • , 
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where 

a n i / \ „9 2 / / s n d , . . d sin 9 d 

A y > := 2K Q (z)rcos9-^ + K (z)rcos9— - k {z) cos 9— + k {z)—^— — , 

9 2 3 9 9 

A (2 ) := 3(K (-z)r cos^) 2 ^-^ + 3«o(#)«o(z)(r cos #) 2 t^; - («o(z) cos ^) 2r g^ + K o( z ) sin<9cos6>— . 

Similarly the two-dimensional Laplacian for any function of (r, 9) is generalized to 

#d#d, 

The perturbation due to the curvature does not affect the two-dimensional Laplacian, and then the 
explicit form of this is calculated as 

d 2 d d 2 



#d#d= — + — + 



dr 2 rdr r 2 d9 2 
2.2 Eigenvalue problems of Laplacians at e = 

The eigenvalue problem associated with the two-dimensional Laplacian, acting on functions of (r, 9) 
that satisfy the Dirichlet condition <3?Ar(a, 9) = 0, is written as 

#d#d$ N {r,9) = -ffir<I> N (r,e). 

An explicit form of the solution is found to be 

* N (r,9) = J n (°^r) e«*, 



and the constant (3n is 

0N = Jn( x q(n)) = 0, 

with N = {n, q(n)},n = 0, ±1, ±2, . . . , and q = 1,2, ... , J n is the n-th order Bessel function and 
x q i n ) is the g-th zero of J n . The orthogonality can be shown as 

¥^$jv#1 =N%$n,m- (4) 

v 

where 

V := { (r,9) | < r < a, < 9 < 2vr }, 
A/% = ITO 2 Jn+l( x q(n))i $N,M '■= 5 n ,m$q,p 

with N = {n, q(n)}, M = {m,p(m)} and S a ,b being the Kronecker delta symbol. 

The eigenvalue problem associated with the three-dimensional Laplacian at e = 0, acting on 
functions of (z,r,9) that satisfy Dirichlet condition (fN^ = at boundary, is solved as 

A^(p N)V (z, r,9) = - 1/3% + (^) 2 | <PN,v(*> r , 0), <Pn,t,(z, r, 9) = $ N {r, 9) sin (J^-z^ . 

Here rj = 1, 2, . . . and the boundary is located at 

z = 0, L, r = a. 

The eigenfunction ipN,ri also satisfies 

(AW + n 2 ) m>1 ,(z, r,9) = - h 2 N + (^) 2 + ^ 2 | <p N , v (z, r, 9), (5) 
with the same boundary conditions where /x is a constant. 
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3 The wave equation in the curved pipe and its naive perturbative 
analysis 

We study the following complex-valued linear partial differential equation for (p(e, t, z, r, 9) € C 

{^^2 -#d#d + A<l>{e,t,z,r,e) = 0, (6) 

where c, [i are constants, the operators appeared here are defined in £j2j and the ranges of coordinates 
are 

-oo < i < oo, < z <L, < r < a, < 6> < 2vr. 
The boundary condition is 

<f>(e,t,z,r,6)\8u = 0, (7) 

where dU denotes spatial boundary. In this paper ([6]) is referred to as the wave equation. 

To obtain an approximate solution we expand the solution based on the lowest order problem 



(e, t, z,r,6)=J2J2 ® N ( r > ^ sin An ^ 



77=1 N 



^( e ,t) = ^(t) + e ^(t) + e ^(t) + o(^), (8) 



and correspondingly 

<f>(e, t, z, r, 6) = (t, z, r, 9) + e</> (1) (t, z, r, 0) + e 2 <^ 2 > (t, z, r, 0) + 0(e 3 ), 
4®(t,z,r,6) := Y,^N(r,e) S m{^z)Af(t), j = 0,l,2,..., 

TV 

where iV := {N,rj} which we call a mode. This expansion of the solution in e and are refereed 
to as the naive perturbation expansion and the mode-amplitude associated with the mode N 
respectively. Note here that this form of the expansion of the solution satisfies the boundary 
conditions, ((TJ), for any value of e. 

To obtain equations for from those for 4>^\ one uses the following identity 




3.1 Unperturbed solution 

To zeroth order in e the equation of motion becomes 

{^|i-A(W}0 (O) (M,r,#)=O, (10) 
where A(°) is defined in ([3]). This can be written in terms of A^ using ([5|) and ([9]) as 
Af + Ofa® =0, := c ] /(^)\^ + ^, 
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where ' denotes the differentiation with respect to t. The solution is 

Af(t) = Af~ ] e~ in ^ + Af +) e in *\ (11) 

where A^.'^ € C are integral constants. 
3.2 First order solution 

To first order in e the equation of motion is written as 

= 2K (z)^ 0) "r cos 6 + K (z)c/> (0)/ r cos 6 - k (z)^^ cos 6 + k (z)^^- 



sm( 



dr 86 

or equivalently, 



2 

Here 



JV"' 



C N,M ■= J *N$Mr cos 0#1 = Tr^n-i + 5 m ,„.+i) jT drr 2 J n (^r) J m 



^(L 2 ) / — f sing d$ M 

C 7VM := / JV ( oTi 1 — C0S6»)#1 



V 



89 r dr 



Km(5 myn -i - 5 m>n+1 ) / drJ n (-^-rj J rn ^-^-r) 
Jo ^ ^ 



•Ep(m) / r i- \ f 1 T 

-7T (d m ,„-l + O m ,n+l) / 

« Jo 



In addition, for a function / of z, we have defined 

(/)s n ,S v := ^ <te sin (^-«) sin Cj^z\ f(z), (f)s v ,C v , ■= dz sin (^r-*) cos (^£z \ f{z). 
Due to (5 n ',n±i in C^N' anc ^ ^n'n' ^ ne right hand side of (fT2l) reduces to 

1(J-) 2 \^fV> A < $(t) = -(—Yv* V V F {1) A {0) (t) 

JV' V 7 »?' n'=n±lq'(n') 

After substituting (fTT]) into (fT2j) one can find the solution (See Sj6]) 

9 / , \ 2 f / A (0 '- ] A {0 ' +) N 

iMaa _ / c I r J d z^ 1 ) + I JV - ' , JV' jn^t 



- \{-£) E E E>^W + 



77' rc'=n±l q'(n') 



(1) A (0 ; -) e -^ + A (0;+) 

+N W F N,n>,q>(n>)tf Q2 _ q2 f> ( 13 ) 

JV JV' J 
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where Rft ft, and Nft ft, are denned as 

R N,N> := N N,N> := 1 ~ fotf.n*/- 

The condition i?^y $7 = 1 at this order is equivalent to the resonance conditions 

Qn,q(n),ri = ^n±l,?'(n±l),r?' • (14) 
On the other hand the condition iV^ = 1 at this order is the case when 

^n,q(ri),ri 7^ ^ra±l,q'(n±l),»/ • 



Note that there are secular terms, oc t, in (|13|) when the resonance condition is satisfied, and that 
the resonance condition, (|14p . does not contain k(z). 

3.3 Second order solution 

To second order in e the equation of motion is written as 

?S- a, M 0,2> 

2K (z)^"rcose + n'JzUW'rcose - k (z)^^- cos9 + k (z) ^- Sm ° 

or oO r 

^2,(0)// , o.. ^NJ f .\/._m2 jlCOV „2/„ N 2n , .2,_A^ (0) 



+3(K O (z)rcos0)> (O)// + 3K O (z)Ko(z)(r cos#)V 0)/ - K^^r-g— cos 2 # + ^(z)-^ cos sin#, 



or equivalently, 



A {2) + O 2 4 (2) 



Here 



H^ eI E £*Mw + E E*M<'>}- < 15 > 

V 7 »/ [n'=n±lg'((i') n'=n,n±2 q> (n') J 



with 



:= / $ N $ Nl r'co S e#l 



v 



* \$n>,n + \^n',n+2 + ^',n- 2 )| dr r 3 J n (^"^ J, 



^(2,2) / ^— 9$AT' 2/)jxi 

C iV JV' := / ^ ^7 r COS 



<9r 



v 

f a f x q{n)\ f x q'(n')\ 

Tri(8 n < >n - 2 ~ 5n',n+2) / drrj n — J J n / — - — J . 



€ i\ 
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One obtains the equation of motion at this order after substituting the lower solutions A^(t) 



and A^(t), given in (jlip and (|13p. into (|15p . The right hand side of (|15p can then be expressed as 



-(— 



E E E ^JV,JV' 

rj' n'=n±l q'(n') 



L V Af N > 



A 



E E E 

rj" n"=n'±lq"{n") 
(0,+) 



F (i) 

N '< N " f 4(0,-) p-in^/,* , 4(0,+) in*,,* 



AT' AT" 



,(2) 



+ L (aAjv) S S ^jV,, 

v 7 7?" n"=n,n±2 ? "(n") 



4(0,-) -if2jv,jv"* _|_ 4(°'+) p^jv.JV"* 



from which one has the solution 



£ \NnJ ^ fa. fa, N > N 'L \Af N , 

vj n'=n±lq'(n') 



E E E 

rj" n"=n'±lq"(n") 



Rfti fin 
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N',N" A (0~) 



A^'TUR 



if 



+ 



JV 



JV 



— iQ^t _|_ AT _ v J. 



N > N " u n* -o 2 

JV JV' 1 



(i) 



_jvv^_ 4 (0,-) I p te e 

" iV iV',JV» 2 _ fl 2 A JV» i ^"^2iJU + ^"ft 2 - ft 2 



JV 



JV' JV 

N '> N " 2in*, A *" i I 40 * + 4ft 2 



JV JV" 



-it 2 



'/V 



JV' 

^'IjV" 4 (0,+) Jo 

+ N '< N " ft 2 ,, -ft 2 ,,, 1 7V ' 7V " 2i% 



JV' JV'' 
2 



JV 

+ iV 



N,N'"-q2 _q2 

Si JV Si JV" 



+ L (aAjv ) ^ ^ ^ 



4(0-) y D 



N - N "ni -o 2 

JV JV" 
£ 6 ^ iv^ 



JV,JV" 



-2ift 



JV 



g— iilj^iit 

'"ft 2 -ft 2 

JV JV'- 



(0,+) 
JV" 



i? 



+ JV 



JV.JV" Q 2 _ r>2 
JV S2 JV" 



(16) 



When the resonance condition at 0(e) is satisfied, the solution at 0(e 2 ) includes the terms being 
proportional to t 2 . From (|16p the resonance conditions between different two modes are obtained 
as 



^n,q(n),r] ~ ^n±2,q" (n±2),rj" 



(17) 



Although when (|17p are not satisfied in addition to the case where (|14p are not satisfied, there are 

(2) 

secular terms in A^ (t) due to the self-mode couplings. The existence of such secular behavior in 



(2) 

A\^ (t) is a notable qualitative difference from the analysis at 0(e). Such secular behavior is due 



v 



S 



to the summation ranges in (|16l) 

£ , and E E = E E • 

n"=n,n±l n'=n±ln"=n'±l n'=n±l n"=n,n±2 

Taking into account the resonance conditions up to 0(e 2 ) one rewrites (|16p explicitly in the 
following cases 

(i) 

n n,q(n), V + ^n±l,g'(n±l),r," for a11 <?'( n ± 1 W and 
n n,q(n), V + ^n±l,g"(n±2),r," ' for a11 ?"( n ± 2 W> 



(ii) 



(iii) 



^n,g(n),r, / ^n±l,g'(n±l),r/' » for all q'(n±l),7]', and 

^n,g(n),r, = ^n±2,g^'(n±2),^' , for some particular <£(n±2),rfi, (18) 



^n,g(n),r, = ^n±l,q{,(n±l),^ ) for some particular gj.(n±l),^, and (19) 

^n, g (n),r, = ^n±2, g ;'(n±2),< , for some particular q"(n±2),T)", (20) 

Here the reason why the case 

±i,g;(n±l),<) for some particular ^(n±l),^, and (21) 

^n,q(n),ri + ^n,q"(n±2), V " , for all q"(n±2),rj", (22) 

does not appear is as follows. If (|2"T|) is satisfied, then one uses ([2Tj) twice, and can find q"(n±2),rj" 
that satisfy 



This is in contradiction to ([22 



^"n,q(n),r) ^"n,q"(n±2),ri' r ' ■ 



3.3.1 Case (i) 

In the case, where both (|14[) and (|17D are not satisfied, one has 

2 



A (o-) A (o,+) 

TV e -iO?ft , jV giH^i 
(1) 



/ v v v f (i) - V i f (2) 



77' n'=n±lg'(n') 



AT' AT 



+ (NS), (23) 



(2) (2 3) 

where F:- T \~ T is real because C / ^ , 6 it vanishes and NS is the abbreviation for the nonres- 

N,N n,q(n),n,q'(n) 

onance terms. The time-dependence of NS is exp(if2^,i) with £1'^, ^Qf] for a given N. 

Thus when the resonance conditions between other different modes are not satisfied one observes 
a secular (divergent in time) behavior, A\J(t) oc t. Correspondingly, our naive perturbative analysis 
could only be valid in a short time range. To improve this naive perturbative result one needs 
another perturbation method that will be discussed in £JH 



9 



3.3.2 Case (ii) 

The solution in the case (ii) is written as 

^ ' rj' n'=n±l g'(n') 

(res) (res) (res) p W ( ,(o-) ,(o,+) 

r," n>'=n>±l q »(n") Q2 N'~ n2 N \~ 2 *% 2i ^N 

x2 (res) (res) (res) / Ao,-) ,(o,+) 

+ Z[jk) ^ S E ^ ^ e " i¥+ 2ifl Le% 1 -(NS), (24) 



where 



(res) (res) (res) 

E E E 

r)" n"=n'±lq"(n") 

are the summations of r]",n" and g" that satisfy (fl~8|) and fi nj9 ( n ) )r; = ^n,<j"(n),?/'> the self-mode 
coupling, for a given set iV. 



3.3.3 Case (iii) 

The solution in the case (iii) is written as 

2 (res) (res) (res) 



tj' n'=n±l q'(n') 

(res) (res) (res) f .(0-) / 2 \ ,(o,+) / 2 \ 

Z_> Z_> ^AT'iV'' "| _OiO.-.. zLO.-. ^ 402. C r 0,0.-.. zLO,. ^ aq2_ 



r," n"=n'±lq"(n") \ N ' \ N NJ N ' \ N 

, X 2(res) (res) (res) / ,(o,-) .(o,+) 

4 {£} E E E ' + ^* I + (- )• c*> 



r/' n"=ri±lq"(n") 



where the summations are performed when (fl~9|) . (f2"0|) and Q n ^ q ( n ) )V = ^ n ,<j"(n),??"; t ne self-mode 
coupling, are satisfied for a given set N. 



4 Amplitude equations 

In this section, we derive amplitude equations by which long time behavior can approximately be 
described for each case in §3.31 

The procedure to obtain an amplitude equation will be shown in the case (i), where there 
is no resonance between different modes. Even in this case as mentioned in §3,3. H there is a 
secular term, oc e 2 t due to the self-mode coupling. This does not make our naive perturbative 
result valid globally in time, and valid only for a short time range. To improve this one uses 
another perturbation method, rather than the naive perturbation method. In this paper we use the 
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renormalization method as a type of systematic perturbation method for removing secular terms 
caused by the use of the naive perturbation method. Using the renormalization method, one can 
obtain the equations of motion that describe the long-time behavior [6] . [7] . [8] . [9] . [10] . In this paper 
we use the method in Ref.pl)] by which reduced equations have been systematically obtained[TT]. 

For the other cases, (ii) and (hi) in £ 13.31 the procedures to obtain amplitude equations are 
similar to the case (i), and we then show only the resultant equations. 



4.1 case (i) 

In the case (i) one concentrates on (I23p . which contains secular terms due to the use of the naive 
perturbation expansion. 

First, to collect secular terms in (123p perturbatively one defines the following polynomials 



N %€ tU N A N ' 

a(0,+) , • 2, n (ren,2) ,(0,+) 



)(ren,2) ,(0,-) 



A 



(pol,+) 
N 



(t) :- 



(26) 
(27) 



where 



N 



m 



N 



C 



N,N' T, 

r\' n'=n±l q'(n') 



(i) 2 

l\n n > ) nl 



F 



(1) 



A 1 



Q2 N,N 
N 



From these definitions, the naive solution can approximately be written as 
Afi{e,t) = Af(t) + eAf(t) + e 2 Af(t) + 0(e 3 ) 



A 



M,-) (t) e -i^t + A f l ' + \t) e 1 ^ 1 + 0(e 3 ) + ( NS ). 



N 



(28) 



Second, one derives the equations which (|26p and (|27p should perturbatively satisfy. From the 
definitions one has 



.4 



(poi,=F) 
N 



{t + r)-A 



N 



(t) 



■ 2 n (ren,2) A (0, T ) 



where r is a real constant. The right hand sides contain A 

terms of A { ? ol ' T \t), not one substitutes A { °' t) 

of the definitions (1261) and (I27D. then 



(o,=F) 



. To obtain the closed equations in 
A^? 01 '^^) + 0{e 2 ) which are the inverse 



.4 



(poZ,=F) 
N 



(t + r)-A^\t) 



Tie 2 n^ n ' 2) Af l > T) (t) + 0(e 4 ). 



N 



Taking the limit, r — > 0, for the both sides one has the renormalization equations 

d_^(ren,^) _ T ,^2n(™. 2 ) ^ (ren,=F) 



dt n 



(29) 



where A^ en '^\t) approximate A^? ol '^\t), and the solutions to (|29|) contain higher order terms in 
e. The explicit forms of the solutions are 



At en ' T \t) = A { p^(0)e^ 2n( r' 2) \ 



N 



(30) 
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Finally one obtains an approximate solution for A^(e, t) that does not contain the secular terms. 
Using j2E]) and ([30]) one has 

The frequency-shift due to the curved boundary is obtained as 



4.2 case (ii) 



We give the renormalization equations for the case (ii). The procedure to obtain these equations 
is same as in §4.11 

To collect the secular terms one defines 



rj n'=n±l(j'(n') 

(res) (res) (res) x p(}) 



sr sr sr n>,n» .( , T ) 

x 2^ Z_> 2^ T 2in*n% Tl -n% * 

r," n"=n'±lq"(n") N N> N 

n / n 2 (res) (res) (res) 

■ r 2 2 f c ^ V V V — F {2) 4 (0 ' T) 

L I jVjv7 2^ 2^ 2^ T 2iQ^ n,n" n" ■ 

\ / ^ n"=n'±l <j"(n") w 

From these equations one can obtain the renormalization equations 

( ren <=F) r. / \ 2 rv / \ 2 



*4r w o 2 



C 



IJV f 2_f/_£_\ \^ \^ \^ F (l) 

77' n' =n±l q'(n') 

(res) (res) (res) ^C 1 ) 

V" V" 1 N',N" Aren,T) 

A 2^ 2^ ^2in«nl-n% x 

r," n"=n'±lq"(n") N N> N 

/ n 2 (res) (res) (res) 

^U£) E E E dsr^4r T) . < 31 > 

V JV 7 rj" n"=n'±lg"(n") N 



Unlike the case (i), the right hand side of (|3ip consists of some different mode amplitudes, then the 

frequency-shift cannot directly be obtained from this form of the renormalization equations. One 
can rewrite (f3Tj) as 

,(ren,=F) _ , • 2 M {H) Arenrf) ( „ 9 \ 

dt A N ~ ±le 2—^ N,N' N' ' W 
N> 

where {M^ 1 ^,} is a real matrix. 
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4.3 case (iii) 

Finally, denning 



/ n 2 (res) (res) (res) pW, 

rf n'=n±l g'(n') 



x 2 (res) (res) (res) n / N 2 

+e 2-f — ] V V V F {1 1 - 
tA \r N ^ ^ N,N> i 

JV/ n' ri=n±Xq'(ri) 



(res) (res) (res) F (}) / 2 



E\ ^ \ ^ -JV'.iV" I ^ t, i (o )=F ) 

v " n"=n'±lq"{n") N \ 

2 (res) (res) (res) 



N 



F ( 2 

rfi n"=n'±lq"(n") ' 



one has 



W4( ren '^) / \2(res) (res) (res) 

tA jy - c ( c ] ^ ^ — Air( 

vf n' '=n±l q> («/) 



dt L \N N J ^ N> 



/ n 2 (res) (res) (res) . . 2 (res) (res) (res) F (}) _ 

I r 2 2 / c \ V V V f C ) V V V AT',iV" 1 ,(ren, T ) 

v iv/ ??' n'=n±l?'(n') V JV 7 r;" n"=n'±l g"(n") ^ JV 
2 (res) (res) (res) 



.2 2 



^ii^j E £ £ -hr i: i,^-- ■ (34) 

L \ N XJ r," n''=n'±lg»(n'') T N 

Here in the course of deriving the renormalization equations we have substituted 

/ x 2 (res) (res) (res) pi}) _ 

v iv 7 r?' n'=n±lg'(n') w 
which are from the definitions of A { ^ ol ' T \t), (33). One can rewrite (134h in the matrix form as 



d 



dt N ^ N > N ' 

N' 



where {M^ n j f/ } is a real matrix, as same as the case of (ii). 



5 Discussion and Conclusions 

In this paper we have explored physical effects due to non-trivial spatial boundary conditions for 
the system of a linear wave equation. As a simple example, we have studied the system with 
Dirichlet type boundary conditions in a prescribed weakly curved pipe. The perturbation scheme 
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was based on the non-perturbed eigenvalue problems of the Laplacians and gave us the set of 
ordinary differential equations. Then it has been observed that self-mode couplings occur and then 
secular terms appear at the second order analysis, in addition to the secular terms due to resonances 
between different modes. The resonance conditions have been derived at each order in the naive 
perturbative analysis and it turns out that these resonance conditions do not contain n(z). Using 
the renormalization method we have derived amplitude equations. In the case where self-mode 
coupling only occurs we have obtained an analytical expression for the frequency-shift due to the 
curved boundary. In the case where there are resonances between different modes, we have derived 
the amplitude equations in a matrix form as renormalization equations. The amplitude equation 
with higher order correction terms can be obtained by applying our procedure straightforwardly 
in each case. Beyond this study, our methodology can be applied to systems with other boundary 
conditions, for example, the prescribed curve is not planer and the cross section is rectangular. In 
the case where the cross section is rectangular, Bessel functions as the eigenfunctions for the two- 
dimensional Laplacian will be replaced with sinusoidal functions. Furthermore, if the given system 
has a weakly nonlinear term, there will be nonlinear resonance terms in the naive perturbation 
expansion in addition to those associated with curved boundary. Thus it is obvious from our 
procedure to obtain the amplitude equations that the corresponding amplitude equations will be 
nonlinear. We believe that our present work and these extensions that follow from this work can 
help elucidate the behavior of systems with non-trivial spacial boundary conditions. 

6 Appendix 

In this appendix we study the following problems 

A + Q 2 A = (a['h + e~ iQt + (a[ +) t + a { +) ) e int , (35) 

and 

A + Q 2 A = {f3[~k + p^) e- in>t + (p[ +) t + ^ +) ) e ify *, (36) 

Here ' denotes the differentiation with respect to t, otj, @j G C and ±f2) € K are given 

constants. 

The solution to (1351) is obtained as 



Similarly, the solution to (|36jl is obtained as 




3 {+) 8 { ~ ] 

w-n' 2 n 2 -n r2 (n 2 -n' 2 ) 
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